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Abstract 



> 

(N 

, Recent observational indications of an accelerating universe enhance the interest in 



studying models with a cosmological constant. We investigate cosmological expan- 
sion (FRW metric) with A > for a general linear equation of state p — wp, w > —1, 



^ I so that the interplay between cosmological vacuum and quintessence is allowed, as 

well. 

bX). Four closed-form solutions (flat universe with any w, and w = 1/3, —1/3, —2/3) 

are given, of which the last one appears to be new. For the open universe a simple 
' relation between solutions with different parameters is established: it turns out that 

5^ , a solution with some w and (properly scaled) A is expressed algebraically via another 

solution with special different values of these parameters. 

The expansion becomes exponential at large times, and the amplitude at the 
exponent depends on the parameters. We study this dependence in detail, deriving 
various representations for the amplitude in terms of integrals and series. The closed- 
form solutions serve as benchmarks, and the solution transformation property noted 
above serves as a useful tool. Among the results obtained, one is that for the open 
universe with relatively small cosmological constant the amplitude is independent 
of the equation of state. Also, estimates of the cosmic age through the observable 
ratio and parameter w are given; when inverted, they provide an estimate 

of w, i. e., the state equation, through the known age of the universe. 
PACS numbers: 98.80.Cq, 04.25.Dm, 04.60-m 
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1 Introduction 



Recent data on the brightness of distant SN la ([1], [2]; see also [3] and the references 
therein), as well as evidence coming from the cosmic age, large scale structure, and cosmic 
microwave background anisotropy combined with the cluster dynamics, indicate, most 
probably, that the observed cosmological expansion is accelerating. Perhaps the most 
natural, although definitely not unique, reason for this acceleration is the presence of a 
cosmic vacuum of nonzero energy and pressure; because of that, investigation of various 
cosmological models including a positive cosmological constant becomes interesting, once 
again. In this paper we study one such model involving matter with an arbitrary (linear) 
equation of state; in particular, an interplay between vacuum and quintessence [4] is 
studied. 

We consider a Friedmann cosmology described by the Friedmann-Robertson- Walker 
metric , 



where dfl is the element of solid angle and k = —1,0, or 1, according to whether the 
universe is open, flat, or closed (we use the units with G = c = 1). The only energy- 
momentum tensor consistent with homogeneity and isotropy of the universe is the one 
corresponding to a perfect fluid, described by the energy density p and pressure p. The 
latter are assumed to satisfy the equation of state 



The range of values of w is chosen on the grounds that w = —1 corresponds to the 
vacuum equation of state, and vacuum contributions are already included in our model 
with the cosmological constant A > 0. Moreover, no physical reasons are known so far 
to justify the values of w which are more negative that the vacuum one, w = —1; it 
is not easy to imagine the kind of physics that would lead to matter with the pressure 
larger than its density. But the following analysis does not depend at all on whether w 
is larger or smaller than one, so we do not limit this parameter from above. Condition 
(|2|) allows for quintessence , that is, for —l<w<0; investigation of its influence on 
cosmological expansion is one of the goals of this paper. 




(1) 



p = wp, 



w = const > — 1 . 



(2) 
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Under these assumptions, Einstein's equations reduce to 

3(l§)' = 8.p + A-34 | = -3^i±i;§, (3) 

\a at J at a at 

the second one expressing energy conservation. By (|2|), this equation is easily integrated 
giving the density in terms of the scale factor, 

where the constant M > characterizes the abundance of matter in the universe. A 
substituion of (^) in the first of equations produces the equation for the scale factor 
a{t) only, 

da\^ 8ttM 1 . A 2 



- 



+ — - A; ; (5) 



dtj 3 a3"'+i 3 

we append it with the Big Bang initial condition 

a(0) = . (6) 

Note that the initial value of the density, p(0), is infinite for w > —1/3, and zero for 
— 1 < u; < —1/3; it is finite and positive for w = —1/3, which value corresponds to the 
so called Einstein's quintessence]^. 

Equation (|5|) with the initial condition (Q) constitute the problem studied below; its 
solution, together with formula (^), specifies completely the corresponding cosmological 
expansion at all times. 

2 Large— Time Expression for the Scale Factor 

2.1 Problem Reformulation and Parameter Range 

For convenience, we rescale the variables according to 

i = V^T, a{t)=[—^\ y{T), (7) 



^The term 'Einstein quintessence' was coined recently in [5] on the grounds of the effective equation 
of state for Einstein's cosmological solution of 1917. 
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and introduce new parameters: 

o / \ \ 2/3(to+l) 

t,^3w + l, u>-2; ^ ^ 4 ( ^ ) = 3 (87rM)-2/('^+2) ^^yu/(.+2) ^ q _ 

A \8nM J ^ ^ ^ ^ 

(8) 

Denoting by a dot the derivative with respect to r, we rewrite (^) and (^) as 

= y-'' + y'^ -ku, y(0) = 0; u > -2; A: = 0,±1; cj > ; (9) 

since we are studying the expansion, a positive square root of the right hand side of (P) 
is, of course, always taken for ij. 

The equation and initial condition are consistent for any value of parameter kuj 
if V is positive, because at small times ~ y^'^ > 0. For = at small times one gets 
~ 1 — kuj, so u is limited by unity for the closed universe {k = 1). When —2 < u < 0, 
the governing equation in the beginning of the expansion becomes ij'^ « —kto, so only 
open and flat universe (A: = — 1, 0) solutions are possible. 

Moreover, we are interested only in the solutions describing an unbounded cosmo- 
logical expansion, i. e., those with ^(t) ^ cxd at r ^ cxd. A simple mechanical analogy 
helps to establish parameter limits ensuring that. Namely, equation (P) describes the 
motion of a particle of a unit mass along the positive y axis, starting at the origin, in 
the potential 

V{y) = -(y-'' + y')<0; 

the total energy of the particle is —kuj. When —2 < u < 0, the potential is monotonically 
decreasing, with the largest value V^(0) = 0. Hence any such motion with a non-negative 
total energy (k = —1,0, open and flat universe) is infinite. The same monotonicity of 
the potential holds for = 0, but the largest value now is V{0) = —1, which yields 
the already known restriction uj < 1 for any infinite motion in the case k = 1 (closed 
universe) . 

Finally, for > the potential is strictly negative and non-monotonic: at y = y^ > 
such that V' (yiy) = it reaches a unique maximum V{yy) = —uj^ given by 

u:, = xmu(y-'' + y^)=y~'' + yl='^yl = ^ \ (10) 
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An infinite motion requires the total energy to be larger than that maximum, —kuj> 
—uJu, which is automatically fulfilled for the open and flat universe; however, for the 
closed universe {k = 1), the parameter lo proves to be limited from above, lo < uJu- This 
constraint has a clear physical meaning: to overwhelm the curvature effect in the closed 
universe and get an unbounded expansion, one has to have enough vacuum repulsion, 
that is, a large enough A. But for > the parameter uj is proportional to a negative 
power of A [see (^], hence a; cannot be too large. 

Note that i/i^ +0 and a;,^ — > 1 + when v — > +0, which allows us to describe the 
cases of positive and zero v uniformly. Also, both Ui, and uji, tend to unity from above 
when ^ oo, hence they have at least one maximum. In fact, each of the functions, 
plotted in Fig.l, has one maximum exactly: max uJi, = 2 at v = 2 (for some reason, 
it is relativistic gas that has the largest 'store of expansion' in the closed universe), 
and maxy^y ~ 1.149 at ~ 7.182. Some values of particular interest are: for u = Q 
{w = —1/3, Einstein's quintessence), yu = 0, = 1; iov u = 1 {w = 0, pressureless 
'dust'), yy = 2-1/3 ^ 0.794, = 3 x 2-^/3 ~ i.ggo; for v = 2 {w = 1/3, relativistic 
gas), yp = 1, uJu = 2; for u = A {w = 1, Zeldovich's superstiff fluid), y^ = 2^/^ w 1.122, 
and again w^, = 3 x 2^2/3 i ggQ_ 

Summarizing, the range of parameters of interest is 

k=l: u > 0, < UJ < uJu] A: = 0,-1: > -2, < < cxd . (11) 

2.2 Expansion Amplitude and Its Basic Properties 

With parameters in the range (0), the cosmological expansion goes indefinitely [unique 
solution to problem (^ grows unboundedly], and at large times it is exponential in time, 
since for y ^ max{l, u} the equation becomes y ^ y [for k = 0, 1 the large time 
condition is apparently relaxed to just y S> 1]. Therefore 

y{T) = y{r,iy,kuj) = A^{kuj)e^[l + o{l)], r ^ oo , (12) 

and we are particularly interested in the behavior of the expansion amplitude A^{kuj) 
for the whole parameter range (|l^). For > and kuj = uj^ (i. e., A: = 1, uj = ujy) the 
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solution y(r) reaches the (unstable) rest point in an infinite time, so Ay{ujy) should 
be zero, which remains to be proved. 

We derive an explicit expression for the amplitude by integrating equation (P), 



T(y) = V, kuj) 



dx 



(13) 



\/ X ^ + x"^ — kuj 

For y, r — > oo the integral here diverges logarithmically, as it should. To get the desired 
expression for the amplitude, one needs to extract this logarithm of y, that is, to add 
and subtract something in the integrand allowing for a closed-form quadrature. The 
most natural choice seems the same integrand with A; = 0, in which case one writes: 



In 



In 



v + 2 



22/(^+2) y 



/-+2)/2 + y^l + y(u+2) 



+ 



\/ X + x"^ — kio y/'. 



X 



dx 



+ 



1 



1 



dx + o(l). 



y ^ oo (14) 



.^J X ^ + x"^ — kuj \J X + 
[unlike (|l^, the integral here converges at infinity]. This expression can be combined 
with the equality ( [T^ ) to write the desired fomula for the amplitude as 

1 1 



\/x + x'^ — kuj \/ X + x2 



dx , 
(15) 

where 

^,(0) = 2-2/(-+2) 

(note that, naturally, 4>v{^) = 0). 

Two features of the expansion amplitude behavior can be established immediately. 
First, when fca; = w — > — 0, the integral (|l5|) for (^yikiS) diverges at x = y^,, so (^uikiS) 
goes to +oo; therefore Ay{ijjy) = 0, as expected. 

Next, we have, for a given z^: 

dlnA^{ktj) d4>y{kijj) 



d{ku)) d{ku)) ' 

and the latter derivative is readily found from ([T5|): 



(16) 



,{kuj) 



d{koj) 



dx 



{x-"" + x"^ - kuj) 



3/2 



> 



(17) 
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Therefore A^{kuj) is a decreasing function of kuj for any given relevant u; in particular, 

A^{-uj) > A^{0)[> A^{u;), u>0,uj<uj^]. (18) 

That means that, as one expects intuitively, for the same parameters w, M, and A, the 
amplitude for the open universe is larger than that for the flat universe which, in its 
turn, is larger than the amplitude for the closed universe if the latter exists. Generally, 
by figuring out the signs of the curvature and parameter and taking into account 
expression (^) for u), one can see that, given k = ±1, u, and the matter abundance M, 
the growth amplitude is larger for larger cosmological constant A. That clearly makes 
physical sense: the larger the contribution of the vacuum repulsion, the stronger the 
expansion is. 

3 Exact Solutions: Flat Universe, Relativistic Gas, 
Einstein Quintessence, Quintessence with w = —2/3, 
and Others 

Exact and, especially, elementary solutions are always most valuable; in the context of 
this study, they specifically shed light on how the expansion amplitude behaves. 

It is not difficult to see that all the integrals involved in the above solution prove to 
be elementary in the following four cases. 

3.1 Flat Universe: k = 



For k = equation (IS) gives the result presented in part in the previous section; here 
is how it looks in full: 



2 

+ 2 



y(-+2)/2 + ^1 + yiu+2) 



sinh ( — — — r 



(19) 



y[T) 

Thus for flat universe the amplitude of the exponential growth of the scale factor is 

^^(0) = 2-2/(^+2)^ 0<A^(0)<1. (20) 

Note also the small time asymptotics expression, 
yir)={^^r) [l + o(l)], 



which is vahd for all the values of kuj and > (the curvature and cosmological constant 
do not play any role in the early universe). For v = and < 0, A; = — 1, the expansion 
becomes linear with the time, in the beginning. 



3.2 Relativistic Gas: w = 1/3 

Here = 2, and we have: 



^ ^ 1 2y2 _ A;cj + 2^2/4 - kuy^ + 1 



2 2-kuj 



< 2 — kuJ , . 2 + kuJ , , kuj 



y{T) = y exp(2r) - exp(-2r) + ^ = V sinh(2r) - ku: sinh^ r , (21) 

so that 

M(k.) = (22) 

(all these results agree with the above ones for k = 0). Note that uj2 = 2 [the definition 
of LUjy is given in ( [To| ) ] , and the expression (^2|) for the amplitude is well defined exactly 
up to a; = 0^2 = 2, for the closed universe {k = 1). 

3.3 Einstein quintessence: w = —1/3 
In this case we have = 0, and 



T{y) = In ^ + — ^+ y ^ yi^^^^ ^ Vl- ku sinh r , (23) 
\'l — kio 



^o(A:w) = , (24) 



hence 



again in agreement with section 3.1. Of course, as seen from (p!0|), cjq = 1. 

3.4 Quintessence with w = —2/3 
Finally, here v = —1, thus 



2V^+1 ' ' ' 

y(r) = ^ ^ ^ exp(r) + ^ [exp(— r) — 1] = V—ku; sinhr + sinh^(T/2) , (25) 
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and 



A-i{kuj) 



1 + 2y/^EJ 



(26) 



Naturally, these expressions remain physically meaningful only for k = —1, 0, coinciding 
with the results of section 3.1 in the latter case of flat universe. Note that solutions 
dl^), (21) and (23) are known (see [6] and the references therein), although perhaps not 
always in exactly the same form; however, solution (p^) appears to be new. 



3.5 Exact Solutions in Terms of Elliptic Functions 

By definition, an elliptic integral is the integral of any rational function of the integration 
variable and the square root of a polynomial of this variable whose degree is three or 
four. Any elliptic integral can be reduced to a combination of elementary functions and 
three standarad elliptic integrals (see c. f. [7]). 



The form of integrals (13) and ( [iq ) determining our cosmological solution and its 
large-time asymptotics prompts an arbitrary power change of the integration variable to 
try reducing them to elliptic integrals. Such an attempt turns out to be successful for 
four different positive values of parameter z^. By the transformation property of section 
5, four reciprocal negative values should be added to the list, making a total of eight 
different exact solutions in terms of elliptic functions. Since these representations are 
rather cumbersome and not very informative, we choose not to give them here explicitly; 
instead, in the table below we show the corresponding values of and w, and the proper 
substitution making the integrals elliptic. 

Elliptic function solutions for > 



V 


w 


Substitution 


6 


5/3 


X = e'^ 


4 


1 


X = r 


1 





X = 


2/3 


-1/9 


x=e 



Note that w = describes pressureless matter ('dust'), and w = 1 corresponds to 
Zeldovich's superstiff fluid [8]. According to (|3^), the reciprocal negative values of u 



9 



allowing for the elliptic function solutions are v = —3/2 — 4/3, —2/3, —1/2, k = —1. 



4 Estimates of the Scale Factor and Expansion Amplitude 

In the general case, no closed-form solution is available, so it is useful to have at least 
some estimates illustrating the behavior of cosmological expansion. Here we derive two 
groups of such estimates differing by the parameters involved in them. It is not difficult 
to check that all the exact elementary solutions satisfy these estimates, as they must. 

4.1 Estimates Independent of the State Equation 

Some estimates of the scale factor and the expansion amplitude follow immediately from 
@. Indeed, we have 



y = yy " + - kuj > yy"^ - kuj 

for all y > when k = —1, 0, and for y > ^fuj when k = 1. From here we obtain: 



^ ^y^-ku?! > 1, y > 0, k = -l; y>V^, k=l. (27) 

In the case of an open universe we integrate this between zero and an arbitrary moment 
of time, which, due to the initial condition (P), results in 



y + V + > exp T . 

Solving this inequality for y provides a simple lower bound of the open universe expansion 
at all times: 

y{T) = y{T,h',—uj) > sinhr, r>0 (resp., y > 0), , k = —1; . (28) 



Sending t ^ oo and comparing to (12), we find the estimate for the growth amplitude: 



Ay{-uj) > V^/2, k = -l . (29) 



In the case of closed universe we integrate inequality (27) up to an arbitrary moment 
of time starting with [see (|l^)] 

JO VX + X'^ — 00 
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with y changing, respectively, from ^/uj to ^(t). Thus we arrive at 



solution of this inequality yields the estimate 

yi^) = y{'^i ^1 ^) > ^ cosh(r - tS), t>t^ (resp., y > ^/uj), k = 1 . (31) 



Again, in the limit r — > cxd this inequality combined with (12) allows for the correspond- 
ing estimate of the expansion amplitude: 



A^{io) > (^/^/2) exp(-r<^), k = I, > . 



(32) 



Note that, by (pO|), t^^ ^ oo when lo io,^ — 0, hence the estimate (|3^) does not 
contradict the amplitude approaching zero at uj^. Clearly, all these estimates do not 
explicitly depend on v, except the dependence on it of r^^ in the case of a closed universe. 



4.2 Estimates Independent of Parameter cu 

On the physical grounds, one expects that the time in which the expansion reaches a 
given value of the scale factor, y, increases from open to flat to closed universe, 

T{y, u, -uj) < T{y, u, 0) [< r(y, u, uj), u > 0] , 

because of the curvature effect. Formally, this is easily derived from the general solution 
([l^ ) by dropping the term ku> in the denominator of the integrand there, and choosing 
the inequality sign accordingly. Using the exact solution ( p!9| ) for the flat universe, we 
thus obtain: 



r < 



1^ + 2 



In 



iv+2)/2 



+ Jl + y(''+^) 



T > 



v + 2 



In 



y(-+2)/2 + ,^1 + y{u+2) 



-1, v>-2- 



k = l, u>0 



(33) 



By inverting these, we get the estimates for the scale factor: 

2/{u+2) 



yir) > 



sinh ( — — — T 



k 



-1, i^>-2; 
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(34) 



sinh ( — — — r 
2 



2/{u+2) 



k = l, 1/ > 



Naturally, the estimates for the growth amplitude implied by ( |34|) are nothing else than 
(ID, in view of (M). 



5 Transformation Property of Cosmological Solutions and 
Expansion Amplitude for the Open Universe 

5.1 Derivation of Transformation Property 

A family of cosmological solutions describing an open universe turns out to have a useful 
transformation property. It relates, in a rather simple algebraic way, a solution with a 
given set of parameters z/, to to another with different values of these parameters. 

To establish this property, let us consider the initial value problem (|9|) with k = —1, 
kuj = —io, and introduce a new unknown function z = z{t) = z{t, u, —uj) instead of 
y = y{T) = yi^^ according to 

y = uj"^ ^2/(.+2) . . ^ ^ ^-./(.+2) . _ (35) 

Since + 2 > 0, the initial condition for z is the same as for 

z(0) = . (36) 
On the other hand, equation (P) after substitution (E^) becomes 



2 ^2 



v + 2 

Introducing the new time. 



f='^r. (37) 



and denoting by a prime the derivative, we rewrite the last equation as 

'\2 



[z) = ^-(-2./(.+2)) ^^2^ ^-{u+2)/2 (38) 
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Clearly, the initial value problem (p8|), ( |3q ) coincides with the original problem (^) 
in which r is replaced with f, v with —2v/{u + 2), and u; with 

^-(.+2)/2^ SO that, by the 

uniqueness theorem, 

z = z{f, V, -Lo) =y[{u + 2)t/2, -2v/{v + 2), -io-^''+''-)l'^^ . 



But then, according to (|35|) and (p 



y(r, u, -co) ^ Lo'/' [y [{u + 2)r/2, -2i./{u + 2), -a;-(-+2)/2^ 



2/{!^+2) 



A: = -1 



(39) 



where the identity sign stresses that the equality holds for all moments of time, as well 



as for any u > —2 and uj > 0. Sending r to infinity provides, by (12), the corresponding 
transformation property of the expansion amplitude, 

2/(i^+2) 



M-U^)^^^' [^-2./(.+2) 



k = -l 



(40) 



In particular, our previous results (23) and ( |2^ ) for = satisfy (|39|) and (|^), respec- 
tively. 

Note also the corresponding transformation of cosmic ages as functions of the scale 
factor y, 



T{y, V, -uj) 



2uj^ 
u + 2 



(i'+2)/2 



-2v/{v + 2), -UJ' 



(41) 



which is easily obtained either by directly inverting (p9|), or by properly changing the 
integration variable in ([l^). 



5.2 Parameter Duality Map 

The correspondence of the indeces involved in 

Mu : ly — > - 



2v 

'u + 2' 



is a one-to-one map of the semiaxis u £ (—2, oo) into itself with the only fixed point 
= 0, so that (0, cxd) is mapped into (—2, 0), and vice versa. Therefore for the case 
of the open universe it is enough to study cosmological solutions with either positive or 
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negative getting then the results for the opposite sign of this parameter immediately 
from ( |39|) and (|40|); we are intensively using this approach below. 

Moreover, it is easy to see that My is a duality map, that is, Ai^ = Ai^^^, or 
My^ =1,1 being the unit map. The whole range of values of u is thus represented as 
a set of pairs of values reciprocal under the map Aiu, with i/ = reciprocal to itself. In 
particular, one reciprocal pair is = 2 (relativistic gas) and = — 1; it is straightforward 
to see that A2{kLo) from (H) and A^i{kuj) from (|6|) satisfy 

Note that for the physical parameter w entering the equation of state, the corre- 
sponding duality relation is 

w + 5/9 

Mw : w — > — — , w> -1 . 

■w + 1 

It gives a one-to-one map of the interval (—1, —1/3) into the semiaxis (—1/3, oo) and 
vice versa; its only fixed point w = —1/3 corresponds to Einstein's quintessence, which 
once again separates one group of cases from the other. The reason for that is well 
known, and it is seen from the equation for the universal acceleration, which is obtained 
by the obvious manipulation with the governing equations (|3|): 

= -^tt{p + 3p) + A = -A7r(3w + l)p + A . (42) 

Therefore matter with w < —1/3 enhances the accelerating effect of vacuum, while 
matter with w > —1/3 diminishes it, and it is only Einstein's quintessence, which is thus 
said to have no gravitational mass, that does not contribute to the acceleration. So, 
given the matter abundance M and k = —1, an expansion with w > —1/3 and certain 
A can in principle be related only to some expansion with w < —1/3 and smaller A, 



as it occurs according to (39). Of course, the Einstein's quintessence solution cannot 



reciprocate with any other, having a different equation of state. 
5.3 Open Universe with a; ^ 1 

As a first example of using the transformation property, let us consider the asymptotic 
behavior of the expansion amplitude for the open universe {k = —1) when uj is large, 
which is easily found from ( pOD and (20). Indeed, in this case a;~('^+^)/^ 0, and we 
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have: 

A^i-u;) = u^'' U^2./(.+2) (0)1 '^^'^'^ [1 + 0(1)] = [1 + 0(1)] . (43) 



This result agrees with ah the exact solutions (|22D, (|24|), and (26), as well as with the 
estimate (p9|). 

Thus the main term giving the large-time behavior of the open universe with u; S> 1 
does not depend on the equation of state. The corresponding full asymptotic series, which 
proves to be converging, is found in section 6.5; all higher order corrections involved in 
it do depend on i. e., on the state equation. 

6 Expansion Amplitude at uj Near the Cutoff Value 

As explained in section 2, cosmological expansion remains infinite, and its growth am- 
plitude is thus well defined, up to u; = ujy for > and closed universe {k = 1), which 
value forms thus a physical cutoff, in this case. For —2 < < and open universe the 
solution exists for all positive values of /cw, so /ccj = should be also treated as a cutoff. 
The expansion amplitude behavior near the cutoff value, which cannot be smooth, is of 
interest in both cases. 

We first consider > 0, = 1, and u) near ojjj, that is, u)y — u)^ +0. In this case the 
main contribution to the integral ( |l5|) for (f>iy{kio) comes from the vicinity of the point 
X = Vv, where 

x'"" + - kuj ^ - UJ + {v + 2){x - yyf . (44) 
Using this approximation it is not difficult to find that 

0i,(6j) = ln(tj;^ - ij) + 0(1) , A; = 1, UJ ^ u)i, - Q , 

and 

A^{uj) = A^{{))K^{uj^-ujfl'^[l + o{l)], k = l, w^w^-O. (45) 

Calculation of the value of the constant Ky requires more tedious analysis carried out 
in section 6.4, which provides 



K,=uj-^'^ei,Y> 



Ln=l 



(46) 
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where the coefficients Cn{i^) = 0(l/n^), n — > oo, are given exphcitly in (|54[). 

For Einstein quintessence, = 0, and any curvature, approximation ( ^4|) works 
uniformly in > only for the function and its first derivative, but not for the second 
one. Therefore, with formally set to zero, it gives a wrong coefficient at the quadratic 
term (2 instead of 1). So for = the right answer is given not by (^) with v = 0, but 
by the exact formula ( p4|) (rederived, by the way, in a different manner in section 7.3). 

Finally, when —2 < z/ < 0, we should consider the open universe only, k = —1, in 
which case we can write, for w — > +0: 



(t>v{kuj) 



1 



1 



X ^ + x"^ + IjJ \/x ^ + 



dx 



UJ 



-l/2u 



di = 0{UJ 



-\l2v. 



[recall that -l/2zy, -(2 + z^)/z^ > 0]. So, by (P) and 



A,(-a;) = A,(0) [i + 0(l^-i/2-)1 , ^,(0) =2-2/(^^+2)^ a; ^ +0 , (47) 



The exact result (26) for A-\(k<jj) agrees with this. Note that the remainder term in (^) 
diverges at w — > +0, hence so does the derivative of Ay{^—Lo\ as expected. Thus, in fact, 
the behavior of the growth amplitude at both cutoffs, uj = lo^, z^ > 0, and liJ = 0, < 0, 
is similar, with the only difference that the limit value of the amplitude is zero in the 
first case and positive in the second. 



7 Power series in ku for Ai,{ku) and its implications 



7.1 Power series in ku for u > 



Repeated differentiation of the equality (^) in kuj provides 
[lnA,(M]^"^ -- 



^W(M, n = l,2,...; /(")(M-§^ 



Also, repeated differentiation of (O) shows that, for z^ > 



4")(0) 



dx 



n Jo {x-" + X 



2\"+l/2' 



n = 1,2,...; 



(48) 



(49) 
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where the Pochhammer symbol (a)„ is defined in a usual way, (a)o = 1, {a)n = 
r{a + n)/T{a) = (q + n — l)(a + n — 2) . . . (q + l)a, n = 1,2, and r(Q) is the 
Euler gamma-function. [Note that the integral in (^) diverges at the lower limit when 
< 0]. Therefore 



In 



A.(0) 



n=l 



nl 



(50) 



and the coefficients of this series are found to be (see Appendix) 

'(0) _ 1 r(a^ + i)r(^ + i; 



n 



1,2,.... 



(51) 



n! 2^/^T n ' 

The radius of convergence and other properties of this series are studied below. Since the 
radius proves to be non-zero, function A^{kuj) is regular at kto = 0. That means that, 
for a given positive z^, the expansion amplitude for the open universe is an analytical 
continuation of that for the closed one, and vice versa. 



7.2 Case u = 2 

The result (|2^ ) for = 2 is easily obtained from the series (|50|). Indeed, by (51) we have 

<^(")(o)_ 1 r(f + i)r(f + 1) _ 1 ^r(n+i) _ i 

n! 2^ nr(n + l) 20rn2"r(n + l) n2"+i' 

where we have used the double argument formula for the gamma- fucntion, 2^''-^r{z)T {z + 1/2) 



Vvrr(2z), with z = (n + l)/2. Thus 
'Aoikuj) 



In 



^o(O) 



1 ^ {koj/2) 
9 2^ 



1 



n=l 



n 



ln(l - kuj/2) , 



and, according to (pO[), it is exactly the formula (p4|): ^2(^1^) = 0.5(2 — /cw)^/^. Note 
that the radius of convergence of the series ( ^0|) coincides with the cutoff value for this 
case, 2 = u;2. 



7.3 Case = 

In this case formula (|5l| ) gives immediately 0i"^(O) = l/2n, so 

' AQ{kuj 



In 



Ao(0) 



1 ^ {kuT 1 1 , ^ 



2 n 

n=l 
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that is, by (pC|), Ao^kuj) = 0.5(1 — kuj)^/'^ , which is exactly the expression (p^. Here 
again the radius of convergence of series (^0|) coincides with the cutoff, 1 = wq, hinting 
that this should probably be a generic property. 



7.4 General case > 



The general result for z/ > comes from studying the behavior of the coefficients of the 
series (pQ) with large numbers n. As demonstrated in the Appendix, 



n — > oo . 



(52) 



This means that the convergence radius of the series (|5^) always coincides, as one would 
now expect, with the cutoff value uju, and allows us to rewrite the series as 



In 



1 



E ^^^^^ + E C„(^) {ku;/u;,r 



n=l 



n 



n=l 



1 



ln(l - kuj/uJ^) + E C'n(^) (^w/^r^)" , 



(53) 



n=l 



where, in view of (51) and (p^). 



1 1 



2 V vr 



2r(^ + i)r(^ + i) 



n! 



n ^ oo 
(54) 



Expression ( |53D is useful for computing Ai,{ku>) in the whole range koj < ujy. Also, for 
A; = 1 it determines completely the asymtotics of the amplitude near the cutoff Wjy, whose 
main term is given above in (^), (^). 



7.5 Series in inverse powers of uj for — 2<z/<0, k = —1 

To obatin some series for Ai^{kuj) with negative (and thus k = —1), we combine the 
above result (^) with the amplitude transformation property (pOD. Writing for brevity 
/i = —2vl(v + 2), the latter provides 



\^A^{-iJ) = lna;^/2_ 



1 



In ^^(-w 



-i//(i/+2)> 



In—: — h- 



2 Vi^T2 



In 
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where we used the value of Afj_{0) from pO|). Since, by definition, fi > 0, we can now 
replace the last term by its series expansion ( |53|) provided that uj^'^/('^+'^^ < uj^^ or 

u + 2 



U X 2/(1^+2) 

2, 



(55) 



[recall [y + 2), {—v~) > 0]. Under this condition we thus obtain the desired series for 
— 2<zv<0, A; = — lin the form [the series coefficients are given in (|5^)]: 



^1/2 1 oo 

lnyl,(-u;) =ln^ + ^=ln[l + K,H] + ^ C7„(^) [^.(a;)]" , (56) 



2u 



2 Vi^T2 



< k^{lo) 



n=l 

-u/2 



< 1 . 



^ + 2' ------ [_^^/(^ + 2)]('^+2)/2 

This series converges if condition (|55| ) holds, and also provides the asymptotic formula 

,,1/2 

^^(-o;) = \l + (-a;('^+2)/2jJ , ^ oo, -2 < 1/ < 0, /fc = -1 , (57) 



which sharpens expression (|43| ) by giving the exact order of the remainder. 



8 Discussion 

We have determined various significant properties of the cosmological solution, particu- 
larly, those of the expansion amplitude, and provided different expressions for the latter. 
The plots of A^{kuj) as a function of ku for several 'popular' values of the parameter i/, 
including positive, zero, and negative ones, are given in Fig. 2. 

An interesting question is what information on the universe can be gotten if the 
value of the expansion amplitude, say, A^,, is known? Even if no other information is 
available, but A^: > 1, from ( pOD and the monotonicity of the amplitude established in 
sec. 2, one immediately concludes that the universe is open, since A,y(kuj) < 1 for k = 1 
(and, of course, u > 0). If, on the other hand, < ^4* < 1, but the equation of state, i. 
e., parameter u, is more or less known, then equation 

A,y{kuj) = A^ 

provides, again due to the monotonic property, a unique solution for ktu, which means 
the sign of the curvature, and a relation between M and A [see (P)]. 
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One way of inferring the amplitude from observational data is as follows. Formula 
(Q) for the density combined with reduces to Svrp = A/y*^"^^, therefore 

Qm _ Stt/o _ 1 

so that 

■ 

This equality holds at all moments of time. If, however, the expansion is in the asymp- 
totic regime, ( |58[ ) and (12) allow for 

Auiku;) « (^^j exp(-r), y(r) » 1 . (59) 

So, if the matter equation of state, the ratio ^a/^m, and the properly scaled [see (0)] 
cosmic age are known, the expansion amplitude is known as well, if the asymptotic stage 
of the expansion is reached. By (|58|), the latter requires 

(^) » ' • <^°> 

which, for not much smaller than negative unity, means just I^a/^^a/ ^ 1- But in 
case the dominant form of matter in the universe is quitessence with the equation of 
state close to the vacuum one, condition ( |60|) might hold, and the expansion be in the 
asymptotic regime, despite the fact that I^a/^^m > 1 is still of the order of unity. 

One can also use (^) and (^) to estimate the cosmic age. If ( |60|) is true [and 
moreover (r^A/^^A/)^^''''^^'' 3> uj for open universe], that is, the large-time regime takes 
place, (|59|) provides 

r ^ l"?r" - ln^^(A:a;) . 

U + Z \Im 

Replacing the expansion amplitude with that of flat universe, Ay{{)), by (|l^) and (pO|) 
leads to 

r<^-ln 4-A , k = -l- -A »max{l, a;}, (61) 



for open, and 
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for closed universe in the asymptotic regime. 

For any moment of the expansion the estimates for the age of the universe come 
from combining expression (^) for the scale factor via observables with 'geometrical' 
inequalities (^): 

r < TFii^A/^M, v), k = -I] T > tfi{^a/0,m, y), k = l . (63) 

Here 



tfi{^a/^m, i^) = — — In 
u + 1 



(64) 



is the exact age of the flat universe with the same parameters. If Q\/Qm, and the 
cosmic age are known accurately enough, its comparison with the bound (|64|) allows for 
an immediate determination of the curvature of the universe. If, instead, only the last 
two parameters and the sign of the curvature are known, the state equation parameter 
can be estimated. 

We now apply the above results to the currently favored cosmological model. It 
includes baryonic and dark matter, both with the same equation of state w = and 
total r^jv/ = 0.3 lb 0.1, and vacuum with i}\ = 0.7 it 0.1. The third component is 
relativistic particles {w = 1/3), but its abundance < 0.001 is so small compared to 
the first two that it should be neglected. This results in exactly the studied model with 
V = \ and r^A/f^A/ = 1-5 -i- 4. Thus 

yshr) =ysh^) =1-1 -1-6-1' (65) 

so, by (|5^), we are still rather far from the large time regime. (However, this is true 
for z/ = 1 or any other v not close to —2; otherwise the scale factor y could be large 
enough. The main reason to assume the pressureless equation of state for dominating 
dark matter is that it is clumped around galaxies, which can happen only if its pressure is 
much smaller than its density. This is clearly true for non-negative pressure, ti; > 0; but 
negative pressure of quintessence with the positive density might provide the clumping 
as well, a question which would be interesting to examine. In case the answer is positive, 
a better determination of the equation of state might become an important observational 
target). 

21 



According to ( |6^ ) and (^), for our universe 

rj^; = 0.7^ 1.0 . (66) 
The dimensional time is related to r by [see (0)] 



T = r W -|- = r (0.7 ^ 0.8) x 10^^ sec = r (21.3 ^ 25.4) Gyr , (67) 

V GrA 

with 

A = Sirpv = 87r(0.7 ± 0.1)pc = (1-0 ^ 1.6) x 10"^^ g/cm^ 

based on 

Pc = (0.7 ± 0.08) X 10^2^ g/cm^ , 

which corresponds to the range of the Hubble constant H = 70 zt W km/ sec ■ Mpc. By 
virtue of (|6^) and (|67|) , for our universe 

Tpi = Tpi (21.3 ^ 25.4) Gyr = (14.9 ^ 17.8) Gyr , (68) 

that is, its age in case the universe is flat. Globular cluster data independently give the 
lower bound of the cosmic age as 

Tu > (12 - 16) Gyr . (69) 

According to these data, our universe is open if Tu = (12 14.9 — 0) Gyr, and it could 
be either open or closed if = (14.9 16) Gyr. Since the first is true not only for 
the centers of both ranges ( p9| ) and (^), but even for the right end of the first and the 
center of the second one, the conclusion is that our universe is most probably open. 

Under this assumption we can now obtain a lower bound (^) of the cosmic age, to 
check the consistency of our estimates. Indeed, using (^), we obtain an estimate of the 
parameter u for our universe as 

.<^=(M5i£)!^!<.,7..„,. (70) 

smh Tu smh r„ 

with Tu ~ 1.54 taken from the middle of the ranges ( |69| ) and ([67|). This and the basic 
expression ( |l3|) for the solution immediately produce 



"(f^A/^M)^/^ dx /-i-i dx 

Vx-'^ + x2 + iVest Jo + 2.7 ^ 
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or, by 

Tu > (9.6 - 11.4) Gyr , 

in nice agreement with (|6S|). 

Finally, to check the consistency of the current cosmological data, we can convert the 
estimate (63) for an open universe into an upper bound for the state equation parameter 



I/, namely: 



ly < 2 I T-^ In 



(71) 



By (6£), ( |67| ) and (^), for our universe this provides 

This is almost entirely consistent with the choosen i/ = 1 (pressureless dust), except for 
a small range of parameters, when the ratio JIa is close to its minimum observed value 
0.6 and simultaneously the cosmic age is close to its maximum of 16 Gyr. 

In conclusion, we point out two possible generalizations of the present study. The first 
of them deals with the phase transition, that is, with a sudden change of the equation 
of state, when w changes abruptly at some moment of time. This case, as well as the 
one with a whole sequence of phase transitions, can be investigated in a fashion similar 
to the above. 

The second way to generalize the study is to consider multicomponent matter, when 

N N 

P=^Pn, P=^Pn, Pn = WnPn, Wn > -I, n = 1, 2, . . . , iV > 1 , 

n=l ra=l 

and the components do not interact. That means that the conservation equation in (^) 
holds for every component separately, yielding pn = Mn/a^^'^"^^\ Mn = const > 0, and 
thus the analog of the governing problem (^) becomes 

N 

f = y-i ^y-^ -kio + Y, Pny-"-, y(o) = o, 

n=2 

with y, T and uj normalized as in (0), (^) using wi and Mi, f„ = 3wn + 1, and 
fin > ^ being the abundance of the nth component normalized appropriately. Writ- 
ing y(T, z^i, z^i, . . . , h'N, kuj) for the solution, we have 

y(r, ku) > y{T, Uj,ku;) 
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for any pertinent j; generally, adding every new component enhances the expansion. 
Assuming vi is the smallest of all the powers, we see that at very large times the expan- 
sion is described by the one-component equation (^) with v = vi; however, the growth 
amplitude depends essentially on the whole expansion history, in other words, on all the 
equations of state. 

As for the other results, all the estimates of sec. 4 remain true, including inequalities 
( p^ ) for the ages of open, flat, and closed universe (of course, the expression for Tpi 
depends now on all the parameters involved). The properly generalized transformation 
property ( |39|) for the open universe solutions also holds; in fact, the number of such 
independent transformations is equal to the number of matter components, N. 
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Appendix. Analysis of Power Series of Ajy{ku] 



First we derive expression (^) for the derivatives (/>i"^(0). By (49), for n = 1,2,... we 
have: 



^ -An) 



(0) 



i/(2n+l) + 2 , 

oo 3,i^(2n+l)/2 I .oo g 2(.+2) ^ 



r{n + 1/2) ' ' Jo (1 + ;r-+2)"+i/2 u + 2 Jo (1 + s)"+i/2 
The last integral is one of the known representations of the beta-function B(x,y) with 
the arguments 

u(2n +l) + 2 un I 1 2n 
x = — — — — — — = — — :t + 77, y = n + - 



2{u + 2) v + 2 2' ^ 2 v + 2' 

which is easily verified by means of the substituion s = t/{\ — t). Therefore 

Y(n + \) ^ + 2 + 2 2'v + 2) 2nv + 2 r(n+i) 
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j_r(_iS + ^)r(^ + i) 

2n 



rfn + i 



and this is, in fact, expression ([50|). 

We now prove the asymtotic formula (^2[) for the coefficients of the seires (^) in the 
general case. To do that, we use the large argument asymptoitcs of the gamma-function, 

z — > +00, 



r(z) = V27r exp[(z- l/2)ln2 - z\ [1 + 0(l/z)] 



and convenient notations 



Z/ + 2' 



Z/ + 2' 



M + = 1 . 



We thus write: 



(0) _ 1 r (nn + 1) r + i) 



n! 2y^n F (n + 1 

and then 

r (un + 1) r (vn + 



r(n + i) 



27r exp(Q„) 



1 + 



n 



n ^ +00 , 



(Al) 



(A.2) 



where 



Q„ = (nn+1/2) ln(nn+l)-(nn+l) + (?jn) ln(?jn+l)-(i;n+l/2)-(n+l/2) ln(n+l) + (n+l) 

(un + 1/2) In(un) + [yn) ln(?;n) — (n + 1/2) Inn + O (1/?^) = 
n ln(n"u^) + (1/2) Inu + O (1/n) . 
Intorducing this into (A. 2), and the resulting equality into (A.l), we end up with 



<^i"^(0) 



1 



^/2^(nV)' 



l + 0(- 

n 



1 {u'^v^'Y 
It remains only to notice that 



l + 0( - 

n 



n ^ +00 . 



(A.S) 



^u + 2j \v + 2j v + 2 

in order to see that (A. 3) is identical to (|52|). 



J 



z/ + 2 ' 
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Figure 1: Parameters for the Closed Universe Solution 
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Figure 2. Ai,{km) for nonnegative v. 
- 1/ = 0, 1 - = 1, 2 - z/ = 2, AS - lower bound 0.5^y{-kLo) 
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Figure 3. A^{km) for negative v: 
l-v = -1/2, 1 - 1/ = -1, 3 - 1/ = -3/2, AS - lower bound {).f>^J{-ku)) 
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